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ABSTRACT 

A new three-dimensional numerical program incorporated with 
comprehensive real gas property models has been developed to sim- 
ulate supersonic reacting flows. The code employs an implicit finite 
volume, Lower-Upper (LU) time-marching method to solve the com- 
plete Navier-Stokes and species equations in a fully-coupled and very 
efficient manner. A chemistry model with nine species and eighteen 
reaction steps are adopted in the program to represent the chemical 
reaction of H2 and air. To demonstrate the capability of the program, 
flow fields of underexpanded hydrogen jets transversely injected into 
supersonic air stream inside the combustors of scramjets are calcu- 
lated. Results clearly depict the flow characteristics, including the 
shock structure, separated flow regions around the injector, and the 
distribution of the combustion products. 
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Hydrogen air mixture passing a 1U ramp. 



ANALYSIS 


• 3 D Navier-Stokes and N a -1 species equations. 

• Brabbs’ finite rate chemistry model. 

- 9 species: H2, H, OH, H2O, O, O2, HO2, H2O2, and N2. 

- 18 reaction steps. 

• Implicit treatment of source terms in species equations. 

• No source term in energy equation for either exothermic or 
endothermic reactions. 

• Temperature and pressure are calculated iteratively from the 
following equations: 
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THERMODYNAMIC AND TRANSPORT PROPERTIES 
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Baldwin-Lomax model. 



NUMERICAL SCHEME 


• Time-marching LU scheme. 

• Implicit treatment of source terms. 


9Q + S- (E* + A n AQ) + -^-(F n + B"AQ) + d 
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NUMERICAL SCHEME 


• Two operators for 3D calculation. 

• Operation count comparable to explicit schemes. 
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PROGRAM VECTORIZATION 


• LHS can be vectorized on planes normal to the sweeping direc- 
tion. 




Lower-Sweep 


Upper-Sweep 
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Resi dual 


CONVERGENCE HISTORY 


• Coarse grid (30x25x30) 



• Present calculation: 

- Residuals reduced by 4 orders of magnitude. 

- Single injection case: (61x39x43) 8 MW, 5 hrs on Cray 2. 

- Dual Injection case: (81x39x43) 11 MW, 6 hrs on Cray 2. 
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1.02.60 x 100 x 50 



MASS-FRACTION CONTOUR ( Two Ports 



53 



FLOW CONFIGURATION 



D=0 . 001 2 






CONCLUDING REMARKS 


b£) 

G 

• H 

+3 

O 

fl 

<D 

Sh 

r O 

CD 

CD 

CX 

C O 

rfl 

bO 


*h 

£ 

03 


£ 

6 W 

G £ 
, O 
<! qG 


U 

0) 

bQ 

• I— ♦ 

>-i 

O 

-M 

o 

0) 

> 

03 

rO 

G w 
G G 
o o 

• rH 

2 "5 

* 2 


CD 


0 

T3 

rfl 

0 

03 

-4^ 

5-i 

a 

CD 

03 

JO 

CO 

+3 

G 

a 

"v 

> 

03 

fl 

fl 

CD 

8 

0 

0 

10 

>> 

CD 

r£> 

-4^ 

pH 

» * 

fl 

fl 

G 

O 

• rH 

CO 

fl 

a 

CO 

CO 

CD 

CD 

CD 

D 

* rH 

CD 

E8 

CD 

<D 

rH 

fl 

>> 

s 

CO 

Jh 

• rH 


CD 

T3 

| 

C 

<D 

> 

CD 

CD 

CO 

CD 

X> 

• rH 

J-H 

rfl 

CO 

<D 

r7-H 


bO 

fl 


CD 

£ 

CO 

CD 

Dfl 

4-^> 

O 

4-^ 


s 

5—1 

o 

G 

03 

a; 
G 
G 

a 

8 JS 

13 -+3 

rG 

H 


G h 

03 £ 
£ G 

+= -+3 

03 G 
rG nz) 

o <3 

* -S 


03 


o 

o 


u 

d 

a s 
6 1 
a 
x 

03 d 
rG 

•+J JjQ 

^ *8 

10 G 

03 
03 
03 

CD o 3 
03 y-l 

3 s 

■»1s 

^ 03 1 
03 G 
+3 


Q 

o 


o 

o 

03 

rG 

-1-3 


bo 

G 


G 

03 

g X 

Jh • H 


O 13 

03 

>> -H 

+3 flj 


bO 

G 


G 

O 


G 

> 

03 

rG 

H 


4-> 

o 
fl 
fl <D 


fl 

O 


u 

a) 

rfl 


f-H 

I 

fl 

o 

fl 


I 

co 

fl 

fl 

Sh 

■4^ 

CO 

-+^ 

CD 

• r— > 

cs 

« 

*4~l 

0 

CO 

1 

cG 

bO 

G 

• H 

o 

fl 

CD 

Sh 


8 
>> G 


G 


03 

(-1 

o m 

*8 

g «* 

1~! +3 

^ £ 

o ^ 

^ o 

'G 'g 

13 G 

• G •-• 

• — i , 

a 10 

CL 03 
ZT -+3 
CO (j 

03 03 

* rH • r— 5 

CD C 

,-. t • rH 

T3 
O 
o 


CD 
CD CO 

*H 

CD 


CO 

fl 

O 

• rH 
+3 

d 

rH 

fl 

o 

*H 

* rH 

o 

CD 

$-4 

<D 

u 

fl 

+3 

o 

fl 

$-4 

4^ 


fl 

o 

• rH 

4^> 

o 

fl 

<D 

*h 

"fl 

O 


8 

. 03 

03 ^ 

O 


r^J 

8 2 

® § 
CD ^ 

fl 

^ .2 

o -+J 

• rH rH 

a _g 

^ £ 
10 t-i 

-13) 

03 

_T_/ • rH 

5-i 10 

$n m 

JV 13 

O -Q 

03 


03 


rG 

H 


G 

03 '"O 

,2 « 

Pn d 

I 


+3 

03 


03 

cH 

03 

rG 

-M 

<4-1 

o 

G 

O 

• rH 

fl 

•4-? 

CD 

fl 

<D 

a 

S-H 

<D 

a 

CD 

CD 

T 3 

co 

G 

O 

• rH 

4-3 

O 

CD 


CD 

O 


^3 

co j 

O 


H 


57 


and more complete combustion of H2 and air. 




